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Abstract 



The production of a stochastic background of relic gravitational waves 
is well known in various works in the literature, where, by using the so 
called adiabatically-amplified zero-point fluctuations process, it has been 
(•**) , shown how the standard inflationary scenario for the early universe can 

■^1- ■ in principle provide a distinctive spectrum of relic gravitational waves. In 

this paper, it is shown that, in general, f(R) theories of gravity produce 
a third massive polarization of gravitational waves and the primordial 
~p^. • production of this polarization is analysed adapting the adiabatically- 

amplified zero-point fluctuations process at this case. In this way, previous 
results, where only particular cases of f(R) theories have been analysed, 
will be generalized. 

The presence of the mass could also have important applications in 
» 1 ■ cosmology, because the fact that gravitational waves can have mass could 

C^ ' give a contribution to the dark matter of the Universe. 

An upper bound for these relic gravitational waves, which arises from 
the WMAP constrains, is also released. 

At the end of the paper, the potential detection of such massive grav- 
itational waves using interferometers like Virgo and LIGO is discussed. 

1 Introduction 

Recently, the data analysis of interferometric gravitational waves (GWs) de- 
tectors has been started (for the current status of GWs interferometers see 
[l][2l|3j|4l[5l[6l[7j|8|) and the scientific community aims in a first direct detec- 
tion of GWs in next years. 



Detectors for GWs will be important for a better knowledge of the Universe 
and also to confirm or to rule out the physical consistency of General Relativity 
or of any other theory of gravitation UHl HU HH H31 HI] ■ In fact, in the context 
of Extended Theories of Gravity, some differences between General Relativity 
and the others theories can be pointed out starting by the linearized theory of 
gravity UM QJEEI]. 

In this picture, detectors for GWs arc in principle sensitive also to a hy- 
pothetical scalar component of gravitational radiation, that appears in ex- 
tended theories of gravity like scalar-tensor gravity and high order theories 
[HI H3 HU HZl HH HH HOI EE IS!], Brans-Dickc theory [53] and string theory 

A possible target of these experiments is the so called stochastic background 
of gravitational waves [25 ] |26 ] 127 ] [28] [29] [30] . 

The production of the primordial part of this stochastic background (relic 
GWs) is well known in the literature starting by the works of [25] [5D] and 
[2ZJ [55]: that, using the so called adiabatically-amplified zero-point fluctuations 
process, have shown in two different ways how the standard inflationary scenario 
for the early universe can in principle provide a distinctive spectrum of relic 
gravitational waves. In [53] [3D] the primordial production has been analysed for 
the scalar component admitted from scalar-tensor gravity. 

In this paper, it is shown that, in general, f(R) theories of gravity produce a 
third massive polarization of gravitational waves and the primordial production 
of this polarization is analysed adapting the adiabatically-amplified zero-point 
fluctuations process at this case. In a recent paper [2T] , such a process has been 
applied to the same class of theories, i.e. the f{R) ones, which will be discussed 
in the present work. But, in |21| a different point of view has been considered. 
In that case, by using a conform analysis, the authors discussed such a process 
in respect to the two standard polarizations which arises from standard General 
Relativity. In the present paper the analysis is focused to the third massive 
polarization. In this way, we generalize previous results where only particular 
cases of f(R) theories have been analysed (see [ID] for example). 

Regarding f(R) theories, even if such theories have to be dismissed either 
because they contradict the current cosmology and solar system tests or generate 
"ghosts" in what would be a quantum version of the finite range gravity |41j. a 
number of authors manage to avoid both of such problems [41] . 

Concerning the presence of the graviton mass, it was considered by many 
authors since Fierz and Pauli (1939) [42] [43]. Important contributions in this 
sense are the ones in [U] 25] HBJ. In these papers, it is emphasized that the 
presence of the mass modifies General Relativity on various scales, which are 
defined by the Compton wavelength of gravitons. This fact leads to dispersion 
of waves and Yukawa-like potentials in the linearized theory. On the other hand, 
it also leads to other modifications in General Relativity solutions in high-order 
(i.e. non-linear) regimes. The fact that gravitational waves can have mass could 
also have important applications in cosmology because such masses could give 
a contribution to the dark matter of the Universe. 

An upper bound for these relic gravitational waves, which arises from the 



WMAP constrains, is also released. About this point, a different interesting 
treatment of the additional polarizations and impact on Cosmic Microwave 
Background has been recently discussed in the good paper [4 7) . 

At the end of the paper the potential detection of such massive gravitational 
waves using interferometers like Virgo and LIGO is discussed. 

2 f(R) theories of gravity 

Let us consider the action 

S= I d 4 x^/^f(R)+£ m , (1) 

where R is the Ricci curvature scalar. 

Equation ([I]) represents the action of the so called f(R) theories of gravity 
[9j [101 HU [HI EH |^1 EH EH ES] in respect to the well known canonical one of 
General Relativity (the Einstein - Hilbert action [HJEZI) which is 

S= f d i xv^gR + £ m . (2) 

The action ((T|) has been analysed in [33J 133 135] in cosmological contexts. 
As we will interact with gravitational waves, i.e. the linearized theory in 
vacuum, C m = will be put and the pure curvature action 

S = Jd i x^f(R) (3) 

will be considered. 



3 Field equations and linearized theory 

Following (32] (note that in this paper we work with 8irG = 1, c = 1 and h = 1), 
the variational principle 

6jd 4 x^df(R)=0 (4) 

in a local Lorentz frame can be used, obtaining: 

f'(R)R^ - \f(R)g^ - /'(iJ) JM! „ + g^af(R) = (5) 

which are the modified Einstein field equations. f'(R) is the derivative of / 
in respect to the Ricci scalar. Writing down, explicitly, the Einstein tensor, eqs. 
(0) become 

<V = jjTgjilaMR) - f'(R)R] + nm-w - g^af(R)}. (6) 



The trace of the field equations © gives 

3D f'(R) + Rf(R) - 2/(12) = 0, (7) 

and, with the identifications [JO] 

*-►/'(*) and w^mnhJimi (8) 

a Klein - Gordon equation for the effective $ scalar field is obtained: 

m 

»-£. (9, 

To study gravitational waves, the linearized theory has to be analysed, with 
a little perturbation of the background, which is assumed given by a near 
Minkowskian background, i.e. a Minkowskian background plus $ = $o (the 
Ricci scalar is assumed constant in the background) [9j [19] . $o is also assumed 
to be a steady minimum for the effective potential V, that called Vq. This as- 
sumption is vital for the further calculations j41j and its physical justification 
arises from the fact that the effective $ scalar field is a function of the Ricci 
curvature and here the linearized theory is developed, i.e. only weak pertur- 
bations near a fixed curvature are considered. Thus, such a minimum has to 
be steady. This is an analysis totally equivalent to the case of the linearization 
process for scalar-tensor gravity, see for example [16]. In Section 4 of [16] the 
authors claim that "we linearize the equations near the background (rj^, ipo ), 
where (po is a minimum of V ". The difference with the present analysis is that 
in [16j the scalar field which is a minimum of the potential arises directly from 
the Brans-Dicke theory [23] , while in the present analysis an effective scalar field 
and an effective potential arise directly from spacetime curvature (the effective 
scalar field is the prime derivative f'(R), see eq. (|8|). 

Thus, the potential presents a square (i.e. parabolic) trend, in function of 
the effective scalar field, near the minimum [41) . i.e. 

V ~ V + )-a8§ 2 =» — ~ m 2 (5$, (10) 

and the constant m has mass dimension. 
Putting 

(11) 

to first order in h^ v and 8$, calling i? Miypcr , R^ and R the linearized quan- 
tities which correspond to R^ upa , R^v an d R, the linearized field equations are 
obtained [HOIEE]: 

Rau - §Vuf = {d^dyhf - TlauOhf) 

(12) 
Dht = m 2 hf, 



where 

Then, from the second of eqs. (|T2|) . the mass can be defined like 



Thus, as the mass is generated by variation of a function of the Ricci scalar, 
in a certain sense, it is generated by variation of spacetimc curvature. In this 
way, the theory is totally generalized for an arbitrary / function of the Ricci 
scalar, improving the results in [9| I19 [ [40 ] . where only particular theories have 
been discussed.. 

R-iivpo and eqs. (fT2j) are invariants for gauge transformations [9l [HI Qjjl [40] 

(15) 
<5$ -> (5$' = 5$; 

then 

^ = hpv - -^r\tiv + Vnvhf (16) 

can be defined, and, considering the transform for the parameter e^ 

□e„ = & 1 V, (17) 

a gauge parallel to the Lorenz one of electromagnetic waves can be chosen: 

9"V = 0. (18) 

In this waythe field equations read like 

□V = ° ( 19 ) 

Ohf = m 2 h f (20) 

Solutions of eqs. (TlT)j) and (|2"U)) are plan waves [T2l[T^] : 

hfj, u = A^ u (jf) exp(ip a x a ) + c.c. (21) 

hf = a(~ft) exp(iq a x a ) + c.c. (22) 



where 



k a = {u^) u=p=\~f\ 



(23) 



q a = (u} m ,J?) u m = yjm 2 + p 2 



In eqs. (fTTJj) and (f^Tj) the equation and the solution for the standard waves 
of General Relativity [31] [32] have been obtained, while eqs. (|20|) and (|22|) 
are respectively the equation and the solution for the massive mode (see also 

HHHmiin]). 

The fact that the dispersion law for the modes of the massive field hf is not 
linear has to be emphasized. The velocity of every "ordinary" (i.e. which arises 
from General Relativity) mode h^ is the light speed c, but the dispersion law 
(the second of eq. (j2"3"]) ) for the modes of hf is that of a massive field which 
can be discussed like a wave-packet HH EH HO]- Also, the group- velocity of a 
wave-packet of hf centred in ~j? is 

n& = ~, (24) 

u 

which is exactly the velocity of a massive particle with mass m and momen- 
tum p. 

From the second of eqs. (f2"3"]) and cq. (IM1) it is simple to obtain: 



\/lo 2 — m 2 . „ x 

u G = ■ (25) 

UJ 



Then, as one wants a constant speed of the wave-packet, it has to be [91 [T2] 



v^ 



(26) 



Now, the analysis can remain in the Lorenz gauge with transformations of the 
type □£„ = 0; this gauge gives a condition of transverse effect for the ordinary 
part of the field: k^A^ = 0, but does not give the transverse effect for the total 
field hav- From cq. (|16p it is 

h^ v = h^,, - -Ti^u + V/fhf. (27) 

At this point, if being in the massless case [HI HH HHJj one puts 

□e'' = 



due»- 



(28) 

5+"/. 

which gives the total transverse effect of the field. But in the massive case 
this is impossible. In fact, by applying the Dalembcrtian operator to the second 
of eqs. (|28p and by using the field equations (fT!?)) and (|20[) it is 

De^ = m 2 h f , (29) 



which is in contrast with the first of eqs. (|28[) . In the same way, it is possible 
to show that it does not exist any linear relation between the tensor field h^ v 
and the massive field hf. Thus, a gauge in which h^ u is purely spatial cannot be 



chosen (i.e. h^o = cannot be put, see eq. (|27|l) . But the traceless condition 
to the field h^ u can be put : 

□e'' = 

(30) 



C M t — 2 . 



These equations imply 



0" V = 0. (31) 

To save the conditions d^h^ and h = transformations like 

□e' 1 = 

(32) 
d^ = 

can be used and, taking ~j? in the z direction, a gauge in which only An, 
A22, and A12 = A21 arc different to zero can be chosen. The condition h = 
gives An = —A 2 2- Now, putting these equations in eq. ([2"T]). it is 

V,(t, z) = A+(t - z)el+J +A x (t- «) e M + / V (t - we*)'?/"'- ( 33 ) 

The term A + (i — z)e/ iv + A x (£ — z)e^ describes the two standard polariza- 
tions of gravitational waves which arise from General Relativity, while the term 
hf(t — vczjrinu is the massive field arising from the high order theory. In other 
words, the function /' of the Ricci scalar generates a third massive polarization 
for gravitational waves which is not present in standard General Relativity. 



4 The primordial production of the third polar- 
ization 

Now, let us consider the primordial physical process, which gave rise to a char- 
acteristic spectrum Q, gw for relic GWs. Such physical process has been analysed 
in different ways: respectively in refs. [551 US] an d [2Z1 US] but only for the 
components of eq. ((33)) which arises from General Relativity. In [29] the pro- 
cess has been extended to scalar-tensor gravity. Actually, the process can be 
further improved showing the primordial production of the third polarization of 
eq. (OH). 

Before starting with the analysis, let us recall that, considering a stochastic 
background of GWs, it can be characterized by a dimensionlcss spectrum [251 

Ha [22 [23 [29] 

n gw (f)^-^, (34) 



where 



3 4 <»> 



is the (actual) critical density energy, p c of the Universe, H the actual value 
of the Hubble expansion rate and dp gw the energy density of relic GWs in the 
frequency range / to f + df. 

The existence of a relic stochastic background of GWs arises from generals 
assumptions, i.e. from a mixing between basic principles of classical theories of 
gravity and of quantum field theory. The strong variations of the gravitational 
field in the early universe amplify the zero-point quantum oscillations and pro- 
duce relic GWs. It is well known that the detection of relic GWs is the only 
way to learn about the evolution of the very early universe, up to the bounds 
of the Planck epoch and the initial singularity [5TJ HS1 HI W\ HI HH] ■ It is very 
important to stress the unavoidable and fundamental character of this mecha- 
nism. The model derives from the inflationary scenario for the early universe 
[35l [36] , which is tuned in a good way with the WMAP data on the Cosmic 
Background Radiation (CBR) (in particular exponential inflation and spectral 
index m 1 [37JI3H])- Inflationary models of the early Universe were analysed in 
the early and middles 1980's (see [35] for a review ), starting from an idea of 
A. Guth [36]. These are cosmological models in which the Universe undergoes 
a brief phase of a very rapid expansion in early times. In this context, the ex- 
pansion could be power-law or exponential in time. Inflationary models provide 
solutions to the horizon and flatness problems and contain a mechanism which 
creates perturbations in all fields. Important for our goals is that this mecha- 
nism also provides a distinctive spectrum of relic GWs. The GWs perturbations 
arise from the uncertainty principle and the spectrum of relic GWs is generated 
from the adiabatically-amplified zero-point fluctuations [21] [25J ES EB ESI EH] . 

Now, the calculation for a simple inflationary model will be shown for the 
third polarization of eq. (|33| , following the works of Allen j25[ [26] that per- 
formed the calculation in the case of standard General Relativity and Corda, 
Capozziello and De Laurentis [29] [30] that extended the process to scalar GWs. 
In a recent paper [21], such a process has been applied to the f(R) theories 
arising from the action ([1]). But, in [5T| a different point of view has been con- 
sidered. In that case, using a conform analysis, the authors discussed such a 
process in respect to the two standard polarizations which arises from standard 
General Relativity. In the following the analysis is focused to the third massive 
polarization. Thus, in a certain sense, one can say that the present analysis is 
an integration of the analysis in [21] . 

It will be assumed that the universe is described by a simple cosmology in 
two stages, an inflationary De Sitter phase and a radiation dominated phase 
[5TJ US] E3 H3 [551 US- The line element of the spacctimc is given by 

ds 2 = a 2 (r))[-dr] 2 + d~t 2 + h^{r 1 ,l?)dx t *dx v ]. (36) 

(|36l) has to be a solution of the general field equations (6). In fact, even if such 
a form is allowed, it could have absolutely different behavior (see [HI H51 B6] V 
For instance one needs to show that inflation is present in the proposed model. 
But, in the case of f(R) theories, which are the ones that we are treating here, 
both of the two conditions are, in general, satisfied. In fact, we recall that the 



original inflation was proposed by Starobinsky in the classical papers [3H1 EH] 
by using the simplest f(R) theory, i.e. the R 2 one. On the other hand, the 
potential presence and the importance of standard De Sitter inflation in the 
general framework of the primordial production of relic gravitational waves has 
been recently shown in |50| considering a different point of view. In that case, 
using a conform analysis, the authors discussed such a process in respect to the 
two standard polarizations which arises from standard General Relativity. In 
the following, the analysis is focused to the third massive polarization. Thus, 
in a certain sense, the present analysis is an integration of the analysis in [SO] ■ 
Remaining in the tapestry of f(R) theories, there are lots of examples in the 
literature where the line element ([55]) , which is the perturbed form of the stan- 
dard conformally flat Robertson- Walker one, results a solution of the general 
field equations (6). One can see for example the recent reviews [5T1 [52l [53] . 

In the line element ()36l) . by considering only the third polarization, the 
metric perturbation (|53")l reduces to 

h^i, = hfl^v, (37) 

where 

10 

V- o 1 o [6S) 

1. 

In the De Sitter phase (77 < 771) the equation of state is P = —p = const, 
the scale factor is 0(77) = r}\r)Q {2r]i — r\)~ x and the Hubble constant is given by 
H(rj) = H ds = crj /r)l. 

In the radiation dominated phase (rj > 771) the equation of state is P = p/3, 
the scale factor is a(rj) = r]/r)o and the Hubble constant is given by H(rf) = 
crjo/n 2 [2T1|251|2S1|22||3JI|. 

Expressing the scale factor in terms of comoving time defined by 

cdt = a{t)dr) (39) 

one gets 

a{t) ex exp(H ds t) (40) 

during the De Sitter phase and 

a(t) oc Vi (41) 

during the radiation dominated phase. The horizon and flatness problems 
are solved if [3SJ 135] 

^ > 10 27 



The third polarization generates weak perturbations h^ u (r], ~af ) of the metric 
(|37l) that can be written, in terms of the conformal time 77, in the form 

V = I^(k)X( V ) exp(t • -£), (42) 

where k is a constant wavevector and 

hf(r},le,-&)=X(ri)expCie--&). (43) 

By putting Y(jj) = a{if)X(rf) one performs the standard linearized calculation 
in which the connections (i.e. the Cristoffel coefficients), the Riemann tensor, 
the Ricci tensor and the Ricci scalar curvature are computed. Then, from the 
Friedman linearized equations, the function Y(rj) satisfies the equation 

Y" + (\t\ 2 -— )Y = (44) 

a 

where ' denotes derivative with respect to the conformal time. Clearly, this 
is the equation for a parametrically disturbed oscillator. 

The solutions of eq. (j44]l give the solutions for the function X(rf) t that can 
be expressed in terms of elementary functions simple cases of half integer Bessel 
or Hankel functions (2TJ [25l [26, 29, 30[ in both of the inflationary and radiation 
dominated eras: 



For 77 < 7/1 

X(n\ = 

0(77) 

for 77 > 7/1 



X(v) = ^rr[l + H ds u- 1 ] exp -ik(r, - m ), (45) 



XWi) = — -—[acxp—ik(r) - 771) + f3 exp ik(n — 771), (46) 

a (v) 

where u) = ck/a is the angular frequency of the wave (that is function of the 
time because of the Constance of k = \ k \ ) , a and j3 are time- independent con- 
stants which can be obtained demanding that both X and dX/ drj are continuous 
at the boundary 77 = 771 between the inflationary and the radiation dominated 
eras of the cosmological expansion. With this constrain it is 



tt=l+i ^m_H^ (47) 

In eqs. (|4"T|) . (|4"5|) u = ck/a(rjo) is the angular frequency that would be 
observed today. Calculations like this are referred in the literature as Bogoliubov 
coefficient methods (2U [25j EH [29l ED] . 
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As inflation damps out any classical or macroscopic perturbations, the min- 
imum allowed level of fluctuations is that required by the uncertainty princi- 
ple. The solution (|4"5")) corresponds precisely to this De Sitter vacuum state 
[HJ |551 HO HE |3D]. Then, if the period of inflation was long enough, the ob- 
servable properties of the Universe today should be the same properties of a 
Universe started in the De Sitter vacuum state. 

In the radiation dominated phase the coefficients of a are the eigenmodes 
which describe particles while the coefficients of (3 are the eigenmodes which de- 
scribe antiparticles. Thus, the number of created particles of angular frequency 
uj in this phase is 

N u = \0 u \ 2 = (^f. (49) 

Now, one can write an expression for the energy spectrum of the relic grav- 
itational waves background in the frequency interval (w, w + duS) as 

d Pgw = 2M^)^ = -j^-- = "i^ry ( 5 °) 

Eq. (|50p can be rewritten in terms of the present day and the De Sitter 
energy-density of the Universe. The Hubble expansion rates is 

tt2 _ 8ttGp c tt2 
n — 3 > n ds 

Then, introducing the Planck density 



PPlanck = T^2~ (51) 



the spectrum is 



n m _ 1 dp sgw _ f dp gw _ 16 p ds 

p c d\af p c df 9 PPlanck 

Some comments are needed. The computation works for a very simplified 
model that does not include the matter dominated era. Including this era, the 
redshift has to be considered. An enlightening computation parallel to the one 
in [25] gives 



n gw {f) = -^^{l + z eq )-\ (53) 



16 pas 

9 PPlanck 

for the waves which at the time in which the Universe was becoming matter 
dominated had a frequency higher than H eq , the Hubble constant at that time. 
This corresponds to frequencies / > (1 + z^y^Ho, where z eq is the redshift 
of the Universe when the matter and radiation energy density were equal. The 
redshift correction in equation (|53[) is needed because the Hubble parameter, 
which is governed by Friedman equations, should be different from the observed 
one Ho for a Universe without matter dominated era. 
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At lower frequencies the spectrum is [HJ [351 HS1 HH1 130] 

JW/)oc/- 2 . (54) 

Moreover, the results (|52p and ([55]) . which are not frequency dependent, 
cannot be applied to all the frequencies. For waves with frequencies less than 
Hq today, the energy density cannot be defined, because the wavelength becomes 
longer than the Hubble radius. In the same way, at high frequencies there is 
a maximum frequency above which the spectrum drops to zero rapidly. In the 
above computation it has been implicitly assumed that the phase transition 
from the inflationary to the radiation dominated epoch is instantaneous. In the 
real Universe this phase transition occurs over some finite time At, and above 
a frequency 

_ o(ti)J_ f55) 

Imax ~ o(t ) At' [ ' 

which is the redshiftcd rate of the transition, VL gw drops rapidly. These two 
cutoffs, at low and high frequencies, to the spectrum force the total energy 
density of the relic gravitational waves to be finite. For GUT energy-scale 
inflation it is [2TJ E3 H3 H3 [3D] ■ 

Pds - 1(T 12 . (56) 



PPlanck 

5 Tuning with WMAP data 

It is well known that WMAP observations put strongly severe restrictions on 
the spectrum of relic gravitational waves. In fig. 1 the spectrum $7 9 ™is mapped 
following [50]: the amplitude is chosen (determined by the ratio — — — ) to be as 
large as possible, consistent with the WMAP constraints on tensor perturbations. 
Nevertheless, because the spectrum falls off oc f~ 2 at low frequencies, this means 
that today, at LIGO- Virgo and LISA frequencies (indicate by the lines in fig. 
1) [20], it is 



^(/)^ioo<9*10- 13 . (57) 

Let us calculate the correspondent strain at s=s lOOif z, where interferometers 
like Virgo and LIGO have a maximum in sensitivity. The well known equation 
for the characteristic amplitude, adapted for the third component of GWs can 
be used [20] : 

h fc (f) ~ 1.26 * lQ-^(±El)^hl 00 n gw (f), (58) 

obtaining [20] 
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h fc (100Hz) < 1.7* 10 -26 . (59) 

Then, as wc expect a sensitivity order of 10 -22 for interferometers at m 
lQQHz, four order of magnitude have to be gained. Let us analyse smaller 
frequencies too. The sensitivity of the Virgo interferometer is of the order of 
10~ 21 at ?s 10Hz and in that case it is [20] 

hf c (WHz) < 1.7*10~ 25 . (60) 

The sensitivity of the LISA interferometer will be of the order of 10 -22 at 
10 -3 sa Hz and in that case it is [2Ti] 

h fc (100Hz)< 1.7* 10 -21 . (61) 

Then, a stochastic background of relic gravitational waves could be, in prin- 
ciple, detected by the LISA interferometer. 

We emphasize the sentence in principle. Actually, one has to take into ac- 
count the Galactic confusion background which will dominate over the instru- 
mental noise. Thus, it is still questionable whether the relic gravitons produced 
in GR could be detected by LISA [41]. 

On the other hand, the assumption that all the tensor perturbations in the 
Universe are due to a stochastic background of GWs is quit strong, but the 
results ([57]) . ([55)1 . (|60[) and ([oT]) can be considered like upper bounds. 

Figure 1, which is adapted from ref. [20j . shows that the spectrum of relic 
SGWs in inflationary models is flat over a wide range of frequencies. The 
horizontal axis is log 10 °f frequency, in Hz. The vertical axis is log 10 Slgsiu- The 
inflationary spectrum rises quickly at low frequencies (wave which re-entered 
in the Hubble sphere after the Universe became matter dominated) and falls 
off above the (appropriately redshifted) frequency scale f m ax associated with 
the fastest characteristic time of the phase transition at the end of inflation. 
The amplitude of the flat region depends only on the energy density during 
the inflationary stage; we have chosen the largest amplitude consistent with the 
WMAP constrains on tensor perturbations. This means that at LIGO and LISA 
frequencies, £l gw {f)h\ m < 9 * 10~ 13 . 

6 Potential detection with interferometers 

Before starting the discussion of the potential interfcrometric detection of the 
massive GWs polarization, there is another point which has to be clarified. 
The attentive reader [5T] asks what would happen to Hulse- Taylor pulsar [S3]. 
As it is clear from eq. ([55]) . by fixing the massless component of the metric 
to GR value, the third term should carry the additional energy which should 
affect evolution of the binary system [?T] . This problem has been discussed by 
Shibata, Nakao and Nakamura in [15] . In such a work, it has been shown that 
the energy associated to the monopole mode of the scalar particle is quite lower 
in respect to the energy carried by ordinary quadrupolc modes arising from 
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Figure 1: The spectrum of relic SGWs in inflationary models is flat over a wide 
range of frequencies. The horizontal axis is log 10 of frequency, in Hz. The 
vertical axis is log 10 Q gS w 
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the linearization of standard General Relativity. Thus, in these conditions, the 
evolution of the binary system should be not affected by the third polarization 
if one remains into observing error bars. 

On the other and, even if the ratio of differential mode to common mode 
is large in the case of the binary pulsar, it could be, in principle, similar and 
potentially small in cosmology. In fact, in the case of f(R) theories such a mode 
arises directly from a function of the Ricci scalar, see the first identification (|SJ), 
i.e. it arises from spacetime curvature, and it is well known that, being the 
production of relic GWs near the bounds of the Planck epoch and the initial 
singularity, in this early Era both of spacetime curvature and its variations were 
very high. Thus, in this cosmological case, the amplitude of the common mode 
could be even dominant in respect to differential modes. 

After these clarifications, let us start the discussion regarding the potential 
detection. 

Even if cqs. (|5^) and (|oT)|) show that the detection of relic GWs is quite 
difficult, people hope in better sensitivity of advanced projects. In this case, 
it is interesting to discuss the interaction between interferometers and massive 
GWs. 

Considering only the third polarization hf(t — VQz)rj^ v the line element as- 
sociated to eq. ([55]) becomes the conformally flat one 

ds 2 = [1 + h f (t, z)]{-dt 2 + dz 2 + dx 2 + dy 2 ). (62) 

As the analysis on the motion of test masses is performed in a laboratory en- 
vironment on Earth, the coordinate system in which the space-time is locally 
flat is typically used and the distance between any two points is given sim- 
ply by the difference in their coordinates in the sense of Newtonian physics 
[T2l IT51 mi [TBI I3T] . This frame is the proper reference frame of a local observer, 
located for example in the position of the beam splitter of an interferometer. 
In this frame GWs manifest themsclf by exerting tidal forces on the masses 
(the mirror and the beam-splitter in the case of an interferometer). A detailed 
analysis of the frame of the local observer is given in ref. J3JJ, sect. 13.6. Here 
only the more important features of this coordinate system are recalled: 

the time coordinate xq is the proper time of the observer O; 

spatial axes are centred in O; 

in the special case of zero acceleration and zero rotation the spatial coor- 
dinates Xj are the proper distances along the axes and the frame of the local 
observer reduces to a local Lorentz frame: in this case the line element reads 



ds 2 = -(dx ) 2 + S ij dx i dx j + 0{\x 3 \ 2 )dx a dx fi . (63) 

The effect of the gravitational wave on test masses is described by the equa- 
tion 



-i&o**. (64) 
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which is the equation for geodesic deviation in this frame. 

Thus, to study the effect of the massive gravitational wave on test masses, 
R l oko have to be computed in the proper reference frame of the local observer. 
But, because the linearized Riemann tensor R^ vpC r is invariant under gauge 
transformations (9j [T2J [13l [14l [31], it can be directly computed from eq. (|37|) . 

From EH it is: 



Rfiuprr — -^{d^dphav + d v d a h^ — dadphfjw — d^dvhap}, 



that, in the case eq. (|37|) . begins 

^o 7 o = ^{d a d h f ?] 0l + dodjhfSQ - d a d 1 hfq m - d d hfS"}; 
the different elements are (only the non zero ones will be written): 

dfhf for a = 7 = 

-d z dthf for a = 3;7 = 



d a d h f r]Q 7 
dod^hfSQ 



d^hj for a = 7 = 



d t d z hf for a = 0;7 = 3 



(65) 



(66) 



(67) 



(68) 



d a d 1 hfr]oo = d a djhf 



-dfh f fc 



a = 7 = 



d\hf for a = 7 = 3 

-d t d z hf for a = 0; 7 = 3 

d z dthf for a — 3; 7 = 

- d d hfS^ = -d 2 z hf for a = 7 . 
Now, putting these results in eq. (|66[) one obtains: 

pi — 1 h 

n aw — 2 / 

"010 = — 2" 1 / 

^030 = \^ h i- 
But, putting the second of equations (|T2|) in the third of eqs. (|7"Tj) . it is 

1 



(69) 



(70) 



(71) 



R\ 



030 



-m ft/, 

which shows that the field is not transversal. 
In fact, using eq. (|64| it results 



(72) 
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hfx, (73) 



2 



V = T^fV (74) 

and 

z = —-m 2 hf{t 1 z)z. (75) 

Then, the effect of the mass is the generation of a longitudinal force (in 
addition to the transverse one). 

For a better understanding of this longitudinal force, let us analyse the effect 
on test masses in the context of the geodesic deviation. 

Following [14] one puts 



-a? 



1_ „o, 1 



R}, j0 = - -d 2 h f (t,z) = --Titfhf + -L ijT n 2 h f . (76) 

m 2 J 

Here the transverse projector with respect to the direction of propagation of 
the GW n, defined by 

Tij = Sij -nSj, (77) 

and the longitudinal projector defined by 

Lij = n&j (78) 

have been used. In this way, the geodesic deviation equation (|64|l can be rewrit- 
ten like 

d 2 1 1 

-j^Xi = -d 2 h f T ljXj - -m 2 h s L ljXr (79) 

Thus, it appears clear that the effect of the mass present in the GW generates 
a longitudinal force proportional to m 2 which is in addition to the transverse 
one. But if v(u) — > 1 in eq. (f25|) one gets m — > 0, and the longitudinal force 
vanishes. Then, it is clear that the longitudinal mode arises from the fact that 
the GW does no propagate at the speed of light. 

Now, let us analyse the detectability of the third polarization computing 
the pattern function of a detector to this massive component. One has to 
recall that it is possible to associate to a detector a detector tensor that, for an 
interferometer with arms along the ue i) directions in respect to the propagating 
gravitational wave (see figure 2), is defined by [H] 

D ij = Uv l v 3 -u l u 3 ). (80) 
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Figure 2: A gravitational wave propagating in the z direction 



If the detector is an interferometer [UI2IH1HG3E1HIIH1III], the signal 
induced by a GW of a generic polarization, here labelled with s(t), is the phase 
shift, which is proportional to [T3] 



s(t) ~ D ij Rio j0 
and, using equations (|76j) . one gets 

s(i) ~ - sin 2 9 cos 2</>. 



(81) 



(82) 

The angular dependence (|82[) . is different from the two well known standard 
ones arising from general relativity which are, respectively (1 + cos 2 6) cos 2</> 
for the + polarization and — cos^su^f? for the x polarization. Thus, in prin- 
ciple, the angular dependence (f82|) could be used to discriminate among f(R) 
theories and general relativity, if present or future detectors will achieve a high 
sensitivity. The third angular dependence is shown in figure. 

For a sake of completeness, let us recall that there is one more problem with 
the potential detection of the common mode. This is because the sensitivity 
curve drawn for different detectors is for the differential mode. The common 
mode is also registered during the GW experiment but it is much noisier, there- 
fore it has different sensitivity level jJT]. Actually, this is correct for studies 
on potential detection in the case of wavelength of the wave much larger than 
the linear dimension of the interferometer (low frequencies approximation). In 
previous analysis, we have implicitly assumed to be in such an approximation, 
and in particular, eqs. (f80|) . (|8Tj) and ([82]) are strictly valid only in this ap- 
proximation. A recent analysis |55j has shown that, in the case in which the 
differential mode is massive, the response function of an interferometer increases 
with increasing frequency differently from the case of massless modes where the 
response function decreases with increasing frequency. This opens important 
perspectives for a potential detection of the massive mode at high frequency. In 



18 



Value 




Theta 



3 



Figure 3: The (dimcnsionless) angular dependence 
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[55] the response function has only been computed in the case of a massive mode 
propagating parallel to one arm of the interferometer, thus, further studies in 
this direction are needed. For example, it will be quite interesting to generalize 
the frequency dependence of the angular pattern (1821 . 

7 Conclusion remarks 

It has been shown that, in general, f(R) theories produce a third massive polar- 
ization of gravitational waves and the primordial production of this polarization 
has been analysed adapting the adiabatically-amplihcd zero-point fluctuations 
process at this case and generalizing previous results in which only particular 
cases have been discussed. 

The presence of the mass could also have important applications in cos- 
mology because the fact that gravitational waves can have mass could give a 
contribution to the dark matter of the Universe. 

An upper bound for these relic gravitational waves, which arises from the 
WMAP constrains, has been also released and at the end of the paper and the 
potential detection of such massive GWs using interferometers like Virgo and 
LIGO has been discussed. 



Acknowledgements 

I would like to thank Francesco Rubanu for helpful advices during my work. I 
have to strongly thank the referees for useful advices and precious comments 
which permitted in improving the paper. 



References 

[1] F. Accrncse ct al. (the Virgo Collaboration) - Class. Quant. Grav. 24, 19, 
S381-S388 (2007) 

[2] C. Corda - Astropart. Phys. 27, No 6, 539-549 (2007); 

[3] C. Corda - Int. J. Mod. Phys. D 16, 9, 1497-1517 (2007) 

[4] B. Willkc et al. - Class. Quant. Grav. 23 8S207-S214 (2006) 

[5] D. Sigg (for the LIGO Scientific Collaboration) 
www.ligo.org/pdf_public/P050036.pdf 

[6] B. Abbott et al. (the LIGO Scientific Collaboration) - Phys. Rev. D 72, 
042002 (2005) 

[7] M. Ando and the TAMA Collaboration - Class. Quant. Grav. 19 7 1615- 
1621 (2002) 



20 



[8] D. Tatsumi, Tsunesada Y and the TAMA Collaboration - Class. Quant. 
Grav. 21 5 S451-S456 (2004) 



[9 
[10 

[11 
[12 

[is; 



[14 
[15 
[16 

[17 
[18 

[19 

[20 
[21 

[22 
[23 
[24 
[25 

[2G 
[27 



C. Corda - J. Cosmol. Astropart. Phys. JCAP04009 (2007) 

C. Corda - Int. Journ. Mod. Phys. A 23, 10, 1521-1535 (2008) 

G. Allemandi, M. Francaviglia, M. L. Ruggicro and A. Tartaglia - Gen. 
Rel. Grav. 37 11 (2005) 

Capozzicllo S and C. Corda - Int. J. Mod. Phys. D 15, 1119 -1150 (2006) 

C. Corda - Response of laser interferometers to scalar gravi- 
tational waves- talk in the Gravitational Waves Data Analy- 
sis Workshop in the General Relativity Trimester of the Insti- 
tut Henri Poincare - Paris 13-17 November 2006, on the web in 
www.luth2.obspm.fr/IHP06/workshops/gwdata/corda.pdf 

C. Corda- Astropart. Phys. 28, 247-250 (2007) 

M. Shibata, K. Nakao and T. Nakamura - Phys. Rev. D 50, 7304 (1994) 

M. Maggiore a nd A. Nicolis - Phys. Rev. D 62 024004 (2000); also in 
lgr-qc/9907055| 



M. E. Tobar , T. Suzuki and K. Kuroda Phys. Rev. D 59 102002 (1999) 

K. Nakao, T. Harada , M. Shibata, S. Kawamura and T. Nakamura - Phys. 
Rev. D 63, 082001 (2001) 

C. Corda and M. F. De Laurentis - Proceedings of the 10th ICATPP Con- 
ference on Astroparticlc, Particle, Space Physics, Detectors and Medical 
Physics - Applications, Villa Olmo, Como, Italy (October 8-12 2007) 

C. Corda - Mod. Phys. Lett. A No. 22, 16, 1167-1173 (2007) 

S. Capozzicllo, C. Corda and M. F. Dc Laurentis - Mod. Phys. Lett. A 22, 
15, 1097-1104 (2007) 

C. Corda - Mod. Phys. Lett. A No. 22, 23, 1727-1735 (2007) 

C. Brans and R. H. Dicke - Phys. Rev. 124, 925 (1961) 

N. Bonasia and M. Gasperini - Phys. Rcw. D 71 104020 (2005) 

B. Allen - Proceedings of the Les Houches School on Astrophysical Sources 
of Gravitational Waves, eds. Jean-Alain Marck and Jean-Pierre Lasota 
(Cambridge University Press, Cambridge, England 1998). 

B. Allen - Phys. Rev. D 37, 2078 (1988) 

L. P. Grishchuk and others - Phys. Usp. 44 1-51 (2001) 



21 



[28] L. P. Grishchuk and others - Usp. Fiz. Nauk 171 3 (2001) 

[29] C. Corda, S. Capozzicllo and M. F. De Laurcntis - AIP Conference Proceed- 
ings, Volume 966, pp 257-263 (2007) - Proceedings of the Fourth Italian- 
Sino Workshop on Relativistic Astrophysics July 20-30 2007, Pescara, Italy 

[30] S. Capozzicllo, C. Corda and M. F. Dc Laurcntis - Mod. Phys. Lett. A 22, 
35, 2647-2655 (2007) 

[31] C. W. Misncr, K. S. Thome and J. A. Wheeler - "Gravitation" - 
W.H.Feeman and Company - 1973 

[32] L. Landau and E. Lifsits - "Tcoria dei campi" - Editori riuniti edition III 
(1999) 

[33] E. Elizalde, S. Nojiri and S.D. Odintsov, Phys. Rev. D70, 043539 (2004) 

[34] G. Cognola, E. Elizalde, S. Nojiri, S.D. Odintsov, L. Scbastiani and S. 
Zerbini, Phys. Rev. D77 (2008) 046009; 

[35] Inagaky T, Nojiri S and Odintsov SD - J. Cosmol. Astropart. Phys. 
JCAP0506010 (2005) 

[35] G. Watson- "An exposition on inflationary cosmology" - North Carolina 
University Press (2000) 

[36] A. Guth - Phys. Rev. 23 347 (1981) 

[37] C. L. Bennett et al - ApJS 148 1 (2003) 

[38] D.N. Spergel et al. - ApJS 148 195 (2003) 

[40] C. Corda Astropart. Phys., 30, 4, 209-215 (2008) 

[41] Private communication with referees 

[42] M. Ficrz and W. Pauli - Proc. R. Soc. A, 173, 211 (1939) 

[43] M. Fierz and W. Pauli - Helv. Phys. Acta 12, 297 (1939) 

[44] A. A. Logunov and M. A. Mcstvirishvili - Thcor. Math. Phys. 65, 971 (1986) 

[45] S. S. Gershtein, A. A. Logunov and M. A. Mestvirishvili - Phys. At. Nucl. 
61, 1420 (1998) 



[46] A. A. Logunov and M. A. Mestvirishvili - |gr-qc/990702T 



[47] D. Bessada and O. Miranda - Class. Quantum Grav. 26 (2009) 045005, also 
in gr-qc/0901.1119 (2009) 

[48] Starobinsky AA - Phys. Lett. B 91, 99 (1980) 

[49] Starobinsky AA - Sov. Phys. JEPT Lett. B 34, 438 (1982) 



22 



[50] S. Capozziello, M. F. De Laurentis and M. Francaviglia - Astropart. Phys. 
2, No 2, 125-129 (2008) 

[51] S. Nojiri and S.D. Odintsov - Int. J. Gcom. Mcth.Mod.Phys.4:115-146 
(2007) 

[52] T.P. Sotiriou and V. Faraoni - larXiv:0805.1"726l 

[53] S. Capozziello and M. Francaviglia - Gen. Rel. Grav. 40, 2-3, (2008) 

[54] R. A. Hulsc and J. H. Taylor, Astrophys. J. Lett. 195, 151 (1975) 

[55] S. Capozziello, C. Corda and M. F. Dc Laurentis - Phys. Lett. B, 669, 5, 
255-259, (2008) 



23 



